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Abstract — An ad hoc network with a finite spatial extent and 
number of nodes or mobiles is analyzed. The mobile locations 
may be drawn from any spatial distribution, and interference- 
avoidance protocols or protection against physical collisions 
among the mobiles may be modeled by placing an exclusion zone 
around each radio. The channel model accounts for the path loss, 
Nakagami fading, and shadowing of each received signal. The 
Nakagami m-parameter can vary among the mobiles, taking any 
positive value for each of the interference signals and any positive 
integer value for the desired signal. The analysis is governed 
by a new exact expression for the outage probability, defined 
to be the probability that the signal-to-interference-and-noise 
ratio (SINR) drops below a threshold, and is conditioned on the 
network geometry and shadowing factors, which have dynamics 
over much slower timescales than the fading. By averaging 
over many network and shadowing realizations, the average 
outage probability and transmission capacity are computed. 
Using the analysis, many aspects of the network performance 
are illuminated. For example, one can determine the influence of 
the choice of spreading factors, the effect of the receiver location 
within the finite network region, and the impact of both the 
fading parameters and the attenuation power laws. 

Index Terms — Ad hoc networks, transmission capacity, Nak- 
agami fading, spread spectrum. 



I. Introduction 

An ad hoc network or peer-to-peer network comprises 
autonomous nodes or mobiles that communicate without a 
centralized control or assistance. Ad hoc networks, which 
have both commercial and military applications, possess no 
supporting infrastructure, fixed or mobile. In addition to being 
essential when a cellular infrastructure is not possible, ad 
hoc networks provide more robustness and flexibility in the 
presence of node or mobile failures than cellular networks. In 
this paper, we evaluate the performance of ad hoc networks of 
independent, identical mobiles that can operate simultaneously 
in the same spectral band. 

In the quest for analytical tractability, most authors (e.g., 
0") - and the many references therein) assume the ad 
hoc network has an infinite number of mobiles spread over 
an infinite area and use the methods of stochastic geometry 
ifTOl . The locations of the mobiles at any time instant are 
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assumed to be spatially distributed as a Poisson point process, 
and performance measures are derived by averaging over the 
spatial distribution. When the network is finite, the Poisson 
point process is an inadequate model primarily because it 
allows an unbounded number of mobiles and does not account 
for the effects of the network boundary. However, this model is 
extremely useful because it enables authors to use Campbell's 
theorem |6|, ITOl to average over the possible network realiza- 
tions and obtain tractable mathematical expressions that vastly 
simplify the performance analysis. Attempts to escape from 
the limitations of the Poisson point process [8| have involved 
restrictive assumptions such as a low density of mobiles ifTTl 
or operation in the high-reliability regime lfl2l . A universal 
assumption is identical fading models for all interference 
and reference signals, for instance Nakagami fading with a 
common Nakagami-m parameter. 

The analysis presented in this paper discards all previous 
assumptions about the spatial distribution of the network 
mobiles and identical fading. In this paper, the spatial extent 
of the network and number of mobiles are finite. Each mobile 
has an arbitrary location distribution with an allowance for the 
mobile's duty factor, shadowing, and possible exclusion zones. 
A new analysis is presented that computes the exact outage 
probabilitjo in the presence of Nakagami fading conditioned 
on the network geometry and shadowing. The Nakagami-m 
parameter can vary among the mobiles, and the reference 
receiver need not be at the center of the network area. 

The main result in this paper is a closed-form expression 
for the conditional outage probability at a reference receiver, 
which is conditioned on the location of the interferers and the 
realization of the shadowing. The expression averages over 
the fading, which has timescales much faster than that of the 
shadowing or node location. The channel from each mobile to 
the reference receiver may have its own distinct Nakagami-m 
parameter, and the ability to vary the Nakagami-m parame- 
ters can be used to model differing line-of-sight conditions 
between the reference receiver and each mobile. The outage 
probability in the presence of Nakagami-m fading has been 
previously considered for interference-limited systems in [13] 
- lfl6l . None of these previous results accommodate thermal 
noise, and none are in closed form because they require the 
computation of derivatives 1131 . Ifl4l or an infinite series lfT31 . 
|[T6l for the most general case of interference with non-integer 
Nakagami-m parameters. Furthermore, these previous works 
assume interference with fixed power and therefore do not 
consider the impact of the spatial model. 



* Outage probability is defined to be the probability that the signal-to- 
interference-and-noise ratio (SINR) drops below a threshold. 
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If unconditional outage probabilities are desired, then the 
conditional outage probability at the reference receiver can be 
averaged over the possible locations of the interfering mobiles 
and over the shadowing realizations. In limited cases, the 
averaging can be done analytically, and the Appendix derives 
the spatially averaged outage probability for the case that 
the reference receiver is at the center of an annular-shaped 
network, the interferers are uniformly distributed, and there 
is no shadowing. The unconditional outage probability can be 
computed for more general cases through the use of Monte 
Carlo simulation, which involves taking the numerical average 
over many network realizations, with the networks drawn 
according to the desired distributions that model the interferer 
locations and the shadowing. Because the fading does not 
need to be realized, such a simulation is extremely fast, with 
execution times that are competitive with the time required to 
compute analytical expressions when they are available. 

A major advantage of our approach over others is that 
any spatial model can be analyzed, and the analysis can 
consider more than mere averages over interfering mobile 
locations. For instance, the probability that a network real- 
ization with specific mobile locations leads to a reference- 
receiver performance that fails to meet an outage constraint 
(the network outage probability) can be determined. Since our 
model allows realistic deployments, it can be incorporated 
into time-dependent random geometric graphs that capture 
the dynamic network status ifPTl . At each sampling instant, 
each mobile receiver has specified parameter values for the 
attenuation power law, the fading, and the shadowing from 
each signal source based on the terrain and mobile locations 
at that instant. The outage probabilities can be calculated by 
applying the closed-form equation presented in this paper. Our 
model alleviates many problems of other analytical models 
such as inability to handle arbitrarily located mobiles, simple 
link models with disc signal coverage, and no modeling of 
border effects. 

The remainder of the paper is as follows. Sectionllllpresents 
a physical interference model of the network. Section [HI] 
presents the derivation and description of the conditional out- 
age probability. Some example outage-probability calculations 
are given in Section [IVJ The issues of spatial averaging and 
shadowing are discussed in Section [VJ while Section [VI] 
examines the effects of mobile density on both the outage 
probability and the transmission capacity. Conclusions are 
drawn in Section IVII1 while the Appendix derives a closed- 
form expression for the spatially averaged outage probability 
in the absence of shadowing when the reference receiver is 
at the center of an annular network of uniformly-distributed 
interferers. 

II. Network Model 

The network comprises M + 2 mobiles that include a 
reference receiver, a desired or reference transmitter Xq, and 
M interfering mobiles X\, ...,Xm- The coordinate system is 
selected such that the reference receiver is at the origin. The 
variable Xi represents both the i th mobile and its location, 
and ||X;|| is the distance from the i th mobile to the receiver. 
The mobiles can be located in any arbitrary two- or three- 
dimensional regions. Two-dimensional coordinates are conve- 



niently represented by allowing Xi to assume a complex value, 
where the real component is the East-West coordinate and the 
imaginary component is the North-South coordinate. 

Although the network area may have any arbitrary shape, 
we assume that it is a circular region of radius r net in the 
subsequent examples. While the reference receiver is always 
located at the origin, the circular region encompassing the 
network does not need to be centered at the origin. Allowing 
the network to be centered at a coordinate other than the origin 
is equivalent to allowing the reference receiver to be located 
away from the center of the network. A circular exclusion zone 
of radius r ex surrounds the reference receiver, and no inter- 
ferers are permitted within the exclusion zone. The exclusion 
zone is a type of guard zone, and a nonzero r ex may be used 
to model the effects of interference-avoidance protocols ifTHl . 
|fl9l , ll20l or a prohibited zone to prevent physical collisions. 
When the reference receiver is located the center of the circular 
network region, the interferers are confined to an annulus with 
inner radius r ex and outer radius r net . For computational and 
analytical convenience, we normalize r net so that r net = 1 
(equivalently, the unit distance is defined to be r ne t). 

Mobile Xi transmits a signal whose average received power 
in the absence of shadowing is Pi at a reference distance do. 
Typically, a common Pi = Pq is used for all mobiles since 
the unknown shadowing conditions make it difficult to adjust 
the power on a per-link basis. Signals may be spread with a 
direct-sequence waveform. Let G denote the processing gain 
or spreading factor of the direct-sequence spread-spectrum 
signal, which is equal to the increase in bandwidth (G = 1 
indicates no spreading). 

In the ad hoc network, the multiple-access communications 
are asynchronous. Since there is no significant advantage to 
using short sequences, long spreading sequences are assumed 
and modeled as random binary sequences with chip duration 
T c . The processing gain G directly reduces the interference 
power. The power of a multiple-access interference signal is 
further reduced by the chip factor h(r ), which is a function of 
the chip waveform and the timing offset r G of the interference 
spreading sequence relative to that of the desired or reference 
signal. Since only timing offsets modulo-T c are relevant, 
< t < T c . In a network of quadriphase direct-sequence 
systems, a multiple-access interference signal with power / 
before despreading is reduced after despreading to the power 
level Ih(T )/G, where ED, E) 

HTo) = ^[RI(To)+R 2 4T c -T )] (1) 

and R^^o) is the partial autocorrelation for the normalized 
chip waveform. Thus, the interference power is effectively 
reduced by the factor G/h(r ). A rectangular chip waveform 
has i?^(r ) = t , and hence 

*..)-l + 2 (*)'-,(*) 

which implies that 1/2 < Ii(t ) < 1. If r D is assumed to 
have a uniform distribution over [0, T c ], then the expected 
value of h(r ) is 2/3. It is assumed henceforth that G/Ii(t ) 
is a constant equal to G/h for all mobiles in the network. 
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After despreading, the power of Xi,i > (i.e., X± is an 
interferer), in the absence of shadowing at distance do is 
Pi = Pi(h/G). For Xq (i.e., the reference transmitter), the 
power after despreading remains Pq = Pq. 

After despreading, Xfs signal is received at the reference 
receiver with an instantaneous power 



p i = P i g i l^l w f{\\X i 



(3) 



where gi is the power gain due to fading, £j is a shadowing 
factor, and /(•) is a path-loss function. The {g,} are indepen- 
dent identically distributed (i.i.d.) and unit-mean with = of, 
where Oj is Nakagami with parameter m t . In Rayleigh fading, 
mi = 1 for all i, and the {gi} are exponentially distributed. In 
the presence of log-normal shadowing, the are i.i.d. zero- 
mean Gaussian with variance a 2 s . In the absence of shadowing, 
£i = 0. For d > do, the path-loss function is expressed as the 
attenuation power law 



f(d) 



(4) 



where a > 2 is the attenuation power-law exponent, and it is 
assumed that do is sufficiently large that the signals are in the 
far field. 

It is assumed that the {g,;} remain fixed for the duration of a 
time interval, but vary independently from interval to interval 
(block fading). With probability p t , the i th interferer transmits 
in the same time interval as the reference signal. The activity 
probability {pi} can be used to model voice-activity factors or 
controlled silence. Although the {pi} need not be the same, it 
is assumed that they are identical in the subsequent examples. 

The instantaneous SINR at the receiver is 



7 = 



A I 



(5) 



m + £ IiPi 



4=1 



where Af is the noise power and 7j is is a Bernoulli variable 
with probability P[I t = 1] = pi and P[l l = 0] = 1 - pi. 
Substituting (O and © into (0 yields 

Po ff olO« o/10 ||X ||- Q 



A I 



J2iiPigiio £i/w \\Xi\\- a 



i=l 

3o^o 



M 



(6) 



r^ 1 +J2 T i9itti 



where T = d^Po/M is the signal-to-noise ratio (SNR) when 
the reference transmitter is at unit distance and fading and 
shadowing are absent, and 

f 10Wi0||x ||- Q i = 

k LtPq 

is the normalized received power due to Xj. 



III. Outage Probability 

Let fl = [fig 5 •••) fijw] represent the set of normalized 
powers given by (0. An outage occurs when the SINR 7 falls 
below an SINR threshold (3, where 7 is given for the particular 
f2 by (|6]l. It follows that the outage probability e for the given 
fj is 

e = P[ 7 <0|n]. (8) 

Because it is conditioned on ft, the outage probability is 
conditioned on the network geometry and shadowing factors, 
which have dynamics over timescales that are much slower 
than the fading. When conditioned on the fading and inter- 
ference, the channel may be assumed to be AWGN, and the 
capacity C(7) of the link conditioned on ft is a function of 
7. The theoretically ideal SINR threshold may be found by 
inverting the capacity function for a given transmission rate 
R, i.e., (3 = C^iR). 

Substituting (0 into ([SJ, and rearranging yields 



e = P 
By defining 



A I 



i=l 

S = f3~ 1 g Qo, Yj = hgiVLi 



A I 



s-E Y < 



1=1 



the outage probability may be expressed as 

e = p [z < r-^n] = f z (r _1 |n) 



(9) 



(10) 



(11) 



(12) 



which is the cumulative distribution function (cdf) of Z con- 
ditioned on fl and evaluated at T _1 . 

Define Fz (z) = 1 — Fz to be the complementary cdf of Z. 
Conditioned on SI, the complementary cdf of Z is 



F z (z\fl) = P[Z > z\fl] = P 



M 

S > z + ^Y^fl 

i=l 



fs(s)ds 



My)dy (13) 



where /y (y ) is the joint probability density function (pdf) of 
the vector (Yi, Ym), the pdf of the gamma-distributed S 
with Nakagami parameter mo is 



fs(s) 



\ fin J 

(mo-iy. 



s mo_1 exp 



(3mos 



(14) 



for s > 0, and the outer integral is over M -dimensional space. 

Successive integration by parts and the assumption that nio 
is a positive integer yield the solution to the inner integral: 



fs(s)ds 



mo — 1 

E 



I3mp 

n 



(15) 
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Defining (3q = /3mo/f2o> ar, d substituting ( fT3T > into (13[ yields 

mo— l 



F 7 (z fl) 



-Poz 



E 

s=0 



xj ... Je-P^^ti + z-^yi) My)dy- (16) 

RM V 1=1 / 

Since s is a positive integer, the binomial theorem indicates 
that 

h-i^§(;>-#) : m 



A multinomial expansion yields 



M 



eh « e (rr 



;>Q \i=l 



0! 



(18) 



where the summation on the right-hand side is over all sets of 
indices that sum to t. Substituting (fTTT i and ( fT8l into ( fl6l ), we 
obtain 



-00* (A>z) 




s=0 
M 



E 

t=0 



Il e ~ PoVl jj )/v(y)dy. (19) 



x E 

Using the fact that the {Y;} are nonnegative and assuming 
that the interfering channels fade independently, 

m o _1 fa \s 
[POZ) 



F z (z\Si) = 



A I 



s! ^ V 

s=0 i=0 
poo A 



'f! 



x E II / 7Te-*»^(»)dy. (20) 

Taking into account the Nakagami fading and the activity 
probability p t , the pdf of fy^y) is 



m 7 - \ 1 



ft| sT 



r(mj) 



f-^e^^ufe) (21) 



where u(y) is the unit-step function, and S(y) is the Dirac 
delta function. Substituting this pdf, the integral in ( l20l i is 

wmc) jy r^ +l j (22) 

where <5^ is the Kronecker delta function, equal to 1 when 
i = 0, and zero otherwise. Substituting (1221 into (|20T > and 
using 

fs\ ftW f si \ ft\\ 1 



(23) 



gives 



mn— 1 s _4 

2 1 



M 

><e n 

ti>Q i=l 



(1 - pi)S ti + 



s=0 t=0 

Pir^ + mt)^ 



1 



(m-i+ii) 



This equation may be written as 

mo — 1 s 



(s-t)\ 



(24) 



(25) 



where 



*i = ( /?o— + 1 ), for % = {1, M}, (26) 



(27) 



it/ 



i t >0 i=l 



Pir(£+rra 
«!r(m») 



for I = 
r i+£ for * > 0. 



(28) 



Equation ( TZTT i may be efficiently computed as follows. 
For each possible t = {0,..,mo — 1}, precompute a ma- 
trix I t whose rows contain all sets of non-negative indices 
{£i, ...,1m} that sum to t. There will be 

<t + m - r 
t 



(29) 



rows and M columns in I t . The 2 t 's may be reused anytime 
the same M is considered. Compute a row vector \l/ containing 
the "J j. For each possible £ = {0, ..,mo — 1}, compute (|281 ) 
using * and place the resulting row into an mo x M matrix G. 
Each term of ( f27T > can be found by using the corresponding 
row from I t as an index into G. Taking the product along 
the length of the resulting row vector gives the corresponding 
term of the summation. More generally, the entire I t matrix 
can be used to index G. To be consistent with matrix-based 
languages, such as Matlab, denote the result of the operation 
as G(l t ). Taking the product along the rows of G(l t ) and 
then the sum down the resulting column gives d27| i. 

IV. Examples 

Consider the network topology shown in the upper right 
corner of Fig. Q] The reference receiver is at the center of 
the network, and M = 28 interferers are within an annular 
region with inner radius r ex = 0.05 and outer radius r net = 1. 
An exclusion zone of radius r ex surrounds each mobile, and 
interfering mobiles are uniformly distributed throughput the 
network area outside of the exclusion zones. Mobiles are 
placed successively according to a uniform clustering model 
as follows. Let X, = r^e^ represent the location of the 
i th mobile. A pair of independent random variables (yi,Zi) 
is selected from the uniform distribution over [0, 1]. From 
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Fig. 1. An example network, which is drawn according to the uniform 
clustering model, is shown in the upper right portion of the figure. The 
reference receiver (indicated by the *) is placed at the origin, and the 
corresponding reference transmitter is located to its right (indicated by the 
• to the right of the receiver). M = 28 interferers (indicated by the •) are 
placed at random within the dashed circle. An exclusion zone surrounds each 
mobile (indicated by solid circles surrounding the mobiles). The figure also 
shows the outage probability as a function of SNR T, conditioned on the 
pictured network topology. Performance is shown for three fading models 
without spreading or shadowing. Analytical expressions are plotted by lines 
while dots represent simulation results (with one million trials per point). 

these variables, the location is initially selected according to 
a uniform spatial distribution over a disk with radius r net by 
setting n = ^/ylr ne t and Qi = 2irzi. If the corresponding Xi 
falls within an exclusion zone of one of the i — 1 previous 
mobile locations, then a new random location is assigned to 
the i th mobile as many times as necessary until it falls outside 
any exclusion zone. 

Once the mobile locations {Xi} are realized, the {Sl^} are 
determined by assuming an attenuation power-law exponent 
a = 3.5, a common transmit power Pi = P for all i, no 
shadowing, and that the reference transmitter is at distance 
H-Xoll =0.1 from the reference receiver. For all results given 
in this paper, the value of pi = 0.5 for all i, and the SINR 
threshold is f3 = dB, which corresponds to the unconstrained 
AWGN capacity limit for a rate-1 channel code and complex 
inputs. 

Example #1. Suppose that spread-spectrum modulation is 
not used (G = h = 1) and that all signals undergo Rayleigh 
fading. Then, m t = 1 for all i, /3 = P/^o = P\\X \\ a , &i = 
\\Xi\\~ a , and d25l l specializes to 

F z (z|n)=e-^n 1 + f°( 1 -^ , (30) 

t— 1 

The outage probability can thus be readily found for any given 
realization of f2. The outage probability, found by evaluating 
d30l > at z = r _1 , is shown along with the spatial locations of 
the mobiles in Fig. Q] Also shown is the outage probability 
generated by simulation, which involves randomly generating 
the mobile locations and the exponentially distributed power 
gains g , ...,<; a/- As can be seen in the figure, the analytical 
and simulation results coincide, which is to be expected 
because ( f30b is exact. Any discrepancy between the curves 
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Fig. 2. Outage probability as a function of SNR T, conditioned on the 
network shown in Fig. [T] with the mixed-fading model. Performance is shown 
both with and without spreading for several values of processing gain (G). 

can be attributed to the finite number of Monte Carlo trials 
(one million trials were executed per SNR point). 

Example #2. Now suppose that the link between the source 
and receiver undergoes Nakagami fading with parameter mo = 
4. The outage probability, found using d25T l through d28l l. is 
also plotted in Fig. [T] The figure shows two choices for the 
Nakagami parameter of the interferers: rrii = 1 and rrii = 4, 
i = l,2,..., M. The rrii — 4 case, denoted by "Nakagami" in 
the figure legend, represents the situation where the reference 
transmitter and interferers are equally visible to the receiver. 
The rrii = 1 case, denoted by "mixed" in the figure legend, 
represents a more typical situation where the interferers are 
not in the line-of-sight. As with the previous example, the 
analytical curves are verified by simulations involving one 
million Monte Carlo trials per SNR point. 

Example #3. The outage probabilities of examples #1 and 
#2 are very high. One way to reduce the outage probability is 
to use spread spectrum. Suppose that direct-sequence spread 
spectrum is used with a processing gain of G and h — 2/3. 
In Fig. [2 the outage probability is shown for direct-sequence 
networks using three different processing gains as well as for 
an unspread network (G = h = 1). A mixed fading model 
(mo = 4 and m, = 1 for i > 1) is used. From this plot, a 
dramatic reduction in outage probability when using direct- 
sequence spreading can be observed. 

V. Spatial Averaging 

Because it is conditioned on the outage probability will 
vary from one network realization to the next. The variability 
in outage probability is illustrated in Fig. [3] which shows the 
outage probability for ten different network realizations. One 
of the networks is the one shown in Fig. Q] while the other 
nine networks were each realized in the same manner; i.e., 
with M = 28 interferers drawn from a uniform clustering 
process with r ex = 0.05 and r net = 1 and the reference 
transmitter placed at distance ||Xo|| = 0.1 from the reference 
receiver. The same set of parameters (a, /?, pi, Pi, and 
used to generate the mixed-fading results of Example #2 were 
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Fig. 3. Outage probability conditioned on several different network 
realizations. A uniform clustering model is assumed with r ne t = 1.0, 
fex = 0.05, M = 28, mixed fading, a = 3.5, and no spreading. Ten 
network realizations were drawn, and the conditional outage probabilities are 
indicated by dashed lines. The average outage probability over 10,000 network 
realizations e is shown by the solid line. 

again used, and signals were not spread. From the figure, it 
can be seen that the outage probabilities of different network 
realizations can vary dramatically. 

In addition to the location of the interferers, fl depends on 
the realization of the shadowing. The shadowing factors {£.;} 
can have any arbitrary distribution and need not be the same 
for all i. In subsequent examples, log-normal shadowing is 
assumed, and the shadow factors are i.i.d. Gaussian with a 
standard deviation a s and zero mean. 

The conditioning on fi can be removed by marginalizing 
Fz(z\ft) with respect to the spatial distribution of the network 
and the distribution of the shadowing. This can be done 
analytically under certain limitations. For instance, with the 
reference receiver at the center of the network, a uniform 
spatial distribution, and no shadowing, the average outage 
probability can be found, as shown in the Appendix. For more 
general cases of interest, such as correlated spatial distributions 
(including uniform clustering processes and Matern processes 
[8)) or the presence of shadowing, closed-form expressions 
are not readily available. However, the outage probability can 
be estimated through Monte Carlo simulation by generating a 
number of f2, computing the outage probability of each, and 
taking the average. Suppose that N networks are generated, 
and let f2 =f2; for the i th network. Then the average outage 
probability may be found from the Monte Carlo estimate of 
the complementary cdf: 

1 N 

Fz(z) = -J2 F z(^)- (3D 

i=i 

The average outage probability is 

e = l-F z (T- 1 ). (32) 

As an example, the solid line in Fig. [3] shows the correspond- 
ing average outage probability averaged over N = 10, 000 
network realizations. 



TABLE I 

Average outage probability when the receiver is at the center 
of the network for various values of m, a, and g with 
uniformly distributed interferers and no shadowing. the 
column marked e s is the average obtained through simulation 
while the column marked e t is the theoretical value. 



M 


a 


G 




e* 


30 
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1 


1.730 x lO" 1 


1.711 x 10" 1 






32 


1.874 x 10 _J 


1.906 x 10 _J 
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1.303 x lO" 1 


1.290 x 10" 1 






32 


3.010 x 10~ J 


3.118 x 10" J 


60 


3 


1 


3.572 x 10- 1 


3.624 x 10" 1 






32 


3.130 x 10~ d 


3.229 x 10" J 




4 


1 


2.518 x 10" 1 


2.574 x 10" 1 






32 


5.429 x 10- J 


5.671 x 



Table J] compares the spatially averaged outage probability 
obtained using the theoretical expression from the Appendix 
against the corresponding value obtained by performing a 
Monte Carlo averaging over N = 10, 000 network realizations. 
For the results shown in Table H] the network radius is 
normalized to r net = 1, the reference transmitter placed at 
distance ||X || = 0.1, the SINR threshold set to /3 = 
dB, and the activity factor set to pi =0.5. An unshadowed 
mixed-fading channel was used with SNR set to T = 10 
dB. An exclusion zone of radius r ex = 0.05 surrounds just 
the reference receiver. Exclusion zones are not placed around 
the interferers or reference transmitter, which ensures that the 
interferers are independently placed according to a uniform 
distribution over the annulus of inner radius r ex and outer 
radius r ne t- Table Q] shows the average outage probability for 
M = {30,60}, a = {3,4}, and G = {1,32}. As can be 
seen, there is a good agreement between the theoretical values 
and those obtained by simulating the interferer locations. The 
Monte Carlo approach has an uncertainty that decreases with 
increasing N, and is an effective approach to handle those 
cases that cannot be solved with the analytical expression of 
the Appendix. 

Although the spatially averaged outage probability e pro- 
vides some insight into the average network behavior, it does 
not capture the dynamics of the spatial distribution. A more 
useful metric for quantifying the spatial variability is the 
probability that the conditional outage probability e is either 
above or below a threshold ex- In particular P[e > er] 
represents the fraction of network realizations that fail to 
meet a minimum required outage probability at the reference 
receiver and can be construed as a network outage probability. 
The complement of the network outage probability P[e < er] 
is the cdf of e: F e {€t)- This probability is shown in Fig. |4] 
for the three fading models without spreading, and for the 
mixed-fading model with spreading. For each curve, the SNR 
is T = 10 dB and the system parameters are the same as 
the corresponding Examples given in Section [IV] Each curve 
was generated by drawing N = 10, 000 networks using the 
same spatial model used to generate the Examples; i.e., with 
M = 28 interferers drawn from a uniform clustering process 
with r ex = 0.05, r net = 1, and no shadowing. For each 
network, the outage probability was computed and compared 
against the threshold et- The curves show the fraction of 
networks whose e does not exceed the threshold. Note that the 
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Fig. 4. Probability that the conditional outage probability e is below the 
threshold outage probability e T in a network with r ne t = 1.0, r ex = 0.05, 
M = 28, and T = 5 dB. TV = 10, 000 network realizations were drawn 
to produce the figure. Results are shown for the three fading models without 
spreading (G = 1) and for the mixed-fading model with spreading (G = 
{8,32, 100}). 




Fig. 5. Probability that the conditional outage probability e is below the 
threshold outage probability e T in a network with r ex = 0.05, r ne t = 1.0, 
M = 28, T = 5 dB, mixed fading, and no spreading (G = 1). N = 10, 000 
network realizations were drawn to produce the figure. Curves shown for no 
shadowing, and shadowing with two values of <J S . 

curves become steeper with increasing G, which shows that 
spreading has the effect of making performance less sensitive 
to the particular network topology. 

Shadowing can be incorporated into the model by simply 
drawing an appropriate set of independent shadowing factors 
for each network realization and using them to compute 
the {Qj} according to (0. In Fig. [5] shadowing was applied 
to the same set of N = 10, 000 networks used to generate 
Fig. 2] Two standard deviations were considered for the log- 
normal shadowing: a s = 2 dB and a s = 8 dB, and again 
a = 3.5. For each shadowed network realization, the outage 
probability was computed for the mixed-fading model without 
spreading (G = 1). All other parameter values are the same 
ones used to produce Fig. |4] The figure shows P[e < ct], 
and the no-shadowing case is again shown for reference. 



TABLE II 

Average outage probability when the receiver is at the center 

(e c ) AND ON THE PERIMETER (e p ) OF THE NETWORK FOR VARIOUS M, a, 
G, AND <J S . 



Parameters 


Outage probabilities 


AI 


Q 


a 


OS 




£p 


30 


3 


i 





0.1528 


0.0608 








8 


0.2102 


0.0940 






32 





0.0017 


0.0012 








8 


0.0112 


0.0085 




4 


1 





0.1113 


0.0459 








8 


0.1410 


0.0636 






32 





0.0028 


0.0017 








8 


0.0123 


0.0089 


60 


3 


1 





0.3395 


0.1328 








8 


0.4102 


0.1892 






32 





0.0030 


0.0017 








8 


0.0163 


0.0107 




4 


1 





0.2333 


0.0954 








8 


0.2769 


0.1247 






32 





0.0052 


0.0027 








8 


0.0184 


0.0117 



As can be seen from the figure, the presence of shadowing 
and increases in a s increase the variability of the conditional 
outage probability, as indicated by the reduced slope of the 
cdf curves. Interestingly, shadowing does not significantly 
alter the average outage probability in this example. For low 
thresholds, such as €t < 0.1, performance is actually better 
with shadowing than without. This behavior is presumably due 
to the fact that, in shadowing, the reference signal power will 
sometimes be much higher than without shadowing. 

In a finite network, the average outage probability depends 
on the location of the reference receiver. Rather than leaving 
the reference receiver at the origin, Table [TT] explores the 
change in performance when the reference receiver moves 
from the center of the radius-r„ et circular network to the 
perimeter of the network. When the reference receiver is at 
the perimeter of the network, the circular network region is 
centered at coordinate —r net . The SNR was set to V = 10 
dB, a mixed-fading channel model was assumed, and other 
parameter values are the same ones used to produce Fig. |5J 
i.e., r ex = 0.05, r net = 1, /3 = dB, and pi = 0.5. The 
interferers were placed according to the uniform clustering 
model and the reference transmitter was placed at distance 
|| Xq II = 0.1 from the reference receiver. For each set of values 
of the parameters G, a, <r s , and M, the outage probability at 
the network center e c and at the network perimeter e p were 
computed by averaging over N = 10, 000 realizations of mo- 
bile placement and shadowing. Two values of each parameter 
were considered: G = {1,32}, a = {3,4}, <r s = {0,8}, 
M = {30, 60}. The table indicates that e p is considerably less 
than l c in the finite network, which cannot be predicted by the 
traditional infinite-network model. This reduction in outage 
probability is more significant for the unspread network, and 
is less pronounced with increasing G. Both e p and e c increase 
as M and a s increase and G decreases. 

As a increases, both e p and e c increase when G = 32, but 
decrease when G = 1. This difference occurs because spread- 
spectrum systems are less susceptible to the near-far problem 
than unspread ones. The increase in a is not enough to cause a 
significant increase in the already high outage probability for 
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Fig. 6. Average outage probability as a function of M for SNR r = 10 dB 
and log-normal shadowing with cr B = 8 dB. Results are shown for the three 
fading models without spreading (G = 1) and for the mixed-fading model 
with spreading (G = {8, 32, 100}). 

unspread systems in those realizations with interferers close 
enough to the reference receiver to cause a near-far problem. 
In the same realizations, the less susceptible spread-spectrum 
systems do experience a significantly increased outage proba- 
bility. 

VI. Number of Interferers and Transmission 
Capacity 

Fig. [6] shows the outage probability averaged over 10, 000 
networks as a function of the number of interferers M in 
the presence of log-normal shadowing with standard-deviation 
cr s = 8 dB. A uniform-clustering network model is again 
assumed with r ex = 0.05 and r net = 1, and the reference 
receiver is at the center of the network (r = 0). The SNR is 
r = 10 dB and the parameters a = 3.5, Pi = 0.5, and (3 = 
dB are again used. Results are shown for the unspread network 
and the three channel models of Examples #1 and #2, as well 
as for the mixed-fading channel model and three spreading 
factors G = {8, 32, 100}. From the curves, it can be observed 
that performance improves with increasing G, and that higher 
values of G reduce the sensitivity to an increased density of 
interfering mobiles. Both of these effects are attributable to 
the improved interference suppression as the processing gain 
increases. 

A trade-off can be seen in Fig. [6] As M increases, the 
performance of each link degrades, yet there will be more 
simultaneous transmissions supported by the network. This 
trade-off can be quantified by the transmission capacity, which 
represents the network throughput per unit area 0~|. It can be 
found by multiplying network density A, which is the number 
of mobiles per unit area, by the rate that bits are successfully 
transmitted over the reference link: 



t = \(l-e)b 



(33) 



where b is the per-link throughput in the absence of an outage, 
in units of bps, and the factor (1 — e) ensures that only 
successful transmissions are counted. 




Fig. 7. Normalized transmission capacity as a function of M for SNR 
r = 10 dB and log-normal shadowing with cr s = 8 dB. Results are shown 
for the three fading models without spreading (G = 1) and for the mixed- 
fading model with spreading (G = {8, 32}). The performance with G = 100 
is not shown and is slightly higher than it is for G = 32. 

Fig. |7] shows the transmission capacity normalized by b 
as a function of M for the same cases shown in Fig. [6] 
As can be seen from the figure, spreading improves the 
transmission capacity due to the corresponding reduction in 
outage probability. The curves become steeper with increasing 
G reflecting the superior ability of spreading to cope with 
increasing network density, provided that enough bandwidth 
is available to support the spread-spectrum signaling. 

VII. Conclusion 

A new and potent method for evaluating the performance 
of an ad hoc network has been presented. Since our model 
allows realistic deployments, it can be incorporated into time- 
dependent random geometric graphs that capture the dynamic 
network status. Factors modeled include a finite spatial extent 
and number of mobiles, any spatial distribution of the mobiles, 
the existence of exclusion zones, path-loss, Nakagami fading, 
and shadowing. The Nakagami m-parameter can vary among 
the mobiles provided that it is an integer for the reference 
transmitter. Using the analysis, many aspects of the network 
performance are illuminated. For example, one can determine 
the influence of the choice of the spreading factor, the effect 
of the receiver location within the finite network region, and 
the impact of the attenuation power law. 

Appendix 

By averaging over Oi = CIm}, the complementary 

cdf of Z is 



(34) 



F z {z) = / fn(u)F z (z\u)du 



where /n(oj) is the joint pdf of fli. Assuming that the {fl;} 
are independent, 



M 



/n(w) =J\fni(ui) 



(35) 
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(^) 



\ mix / 



dx 



(39) 



m -l 



A/ 



F z (z) 



E (/M S E 



s=0 



t=0 



(s t) 



e n 



d>0 i- 



(1 - Pi )S ti + 



2p i T(£ i + mi )m?> [J ( Ci r% et ) - J (*t%)] 



\rl et -rl x ){imrni)Pt i+li) (™i 



-) 



(44) 



where /Oj(wi) is the pdf of Oj. Define 15^ = ,« > 1, to 
be the inverse normalized power of the i th interferer. Using 
d35l l and the definition of I3i, ( f34b may be expressed as 



Fz(z) = 



/ (jijuM^j F z (z\{x^\ ..x^})dx. (36) 



If the location of interferer Xi is uniform within the annulus 
of inner radius r ex and outer radius r net , then its distance 
r'i = \\Xi\ \ from the origin has pdf 



2r 



net ' ex 



o 



for r e . T < r < r net 
elsewhere. 



(37) 



From the definition of I5 l and ©, = (G/h)(P Q /Pi)rf,i > 
1, in the absence of shadowing, and it follows that the pdf of 

Vi is 



2x 



elsewhere 



for ar% x <x< ar% et 



(38) 

where a = (G/h)(P /Pi). Substituting d25j and d38j into 
(|36"1 > and extracting factors from the integral that are indepen- 
dent of the variable of integration yields equation ( |39l at the 
top of the page, where /3q = f3mo\\Xo\\ a in the absence of 
fading. 

By factoring the denominator and splitting the integration 
interval, the integral in d39l is 



b X 



(^) 



\ TYliX J 



V ' 1 dx = -4^-^ [I (b) - I (a)} (40) 



(m,+<.) 



5 



where a = Cir" x , b = Cir" et , and 



I(y)= / x 
Jo 



v / \ -(mi-MO 

(m.+l-i) 



dx. (41) 



By performing the change of variable x = yv, 



I(y)= / (j/ i y)( m<+ «- 1 ) (1 + 



o 



uiiyv 



Po 



ydv 



V 

7 2T2^1 

(m,: + 4 ) 



H hi; — — 



(42) 



where 2-P1 is the Gauss hypergeometric function, which has 
the integral representation: 

2 F 1 ([a,b];c]x) = 

r( c ) 



r(6)r( c -6) Jo 



.6-1/ 



(1 - v) c -"- L (l -xv)~ a dv 
(43) 



and c = b + 1, T(c)/(T(b)T(c — b)) = b. The hypergeometric 
function can be represented by an infinite series if \x\ < 1 and 
is evaluated by analytical continuation if \x\ > 1 ||231 . l24l . 
The function is widely known and is implemented as a single 
function call in most mathematical programming languages, 
including Matlab. The complementary cdf of Z can be found 
by substituting (l40l > into (|39l , with the /(•) function as defined 
in (l42l . The result is equation ( l44b at the top of the page, where 



J(y) = 2F1 



m l +l l7 m i + - 



2 
o 



Po 



(45) 
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